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Simulations that use the clover action in quenched QCD calculations have a lower limit to the quark mass that 
can be reached, because of the fluctuations caused by exceptional configurations. From this low statistics study, 
we find that the twisted clover action, recently introduced by the ALPHA collaboration, can be used to simulate 
quenched QCD at quark masses below those attainable by simulations that use the clover action. 



1. INTRODUCTION 

The masses of the so called light quarks used in 
lattice QCD calculations are not light enough for 
chiral perturbation theory to be a reliable guide 
to extrapolating the results to the physical quark 
masses. For example, in our recent work on de- 
termining the parameters of the chiral QCD La- 
grangian Q , we were unable to test fully the chi- 
ral perturbation theory predictions because the 
dynamical quarks were too heavy and there was 
a lower limit on the mass of the valence quarks 
caused by exceptional configurations. 

The ALPHA collaboration have developed || 
the twisted clover fermion action, specifically to 
study the light quark mass region of QCD. The 
ALPHA collaboration || have used the twisted 
clover action in numerical simulations within the 
Schrodinger functional formalism. We report the 
first computation of hadron spectroscopy using 
the twisted clover fermion operator with periodic 
in space and anti-periodic in time boundary con- 
ditions. 

2. THE TWISTED-FERMION OPERA- 
TOR 

S F = ^2ip{x)(D + m q + in q j 5 T 3 )ip(x) (1) 



The operator D is the fermion operator that is 
independent of the quark mass. The mass term 
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(2) 



has no additive renormalisation, hence it provides 
a lower limit to the eigenvalue spectrum of the 
fermion operator. 

The main disadvantages of the twisted clover 
action are that parity and flavour symmetry are 
explicitly broken. However, the flavour symmetry 
breaking looks very different to the mechanism 
in the staggered fermion formulation so the 
effects of the flavour symmetry breaking may be 
less of a problem. 

The transformation 



y/ -> e <a75T3/2 yj 



ijj' — > tpe 



«*75l"3/ 2 



(3) 



transforms the action in Eq. [I], when D anti- 
commutes with 75 as in the continuum, into the 
standard QCD action 



S F =Y,i>'{*){D + M Q )i)'{x) 

X 

where 



(4) 



Mn = 



tana 



m„ 



*Presented by C. McNeile. 



The ALPHA collaboration ^,[| advocate the fol- 
lowing procedure. The simulations are carried 
out using the twisted clover operator on the lat- 
tice. The results are then matched to the contin- 
uum twisted fermion operator. The rotations in 
Eq. H are then used to rotate back to the standard 
QCD fermion action in the continuum. The rota- 
tions should reduce the effect of the flavour and 
parity symmetry breaking on the physical spec- 
trum. The connection between QCD and twisted 
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QCD could be more complicated, when there 
is spontaneous breakdown of parity and flavour 
symmetry, such as in the Aoki phase ||. 

Another way to look at the connection between 
the twisted and the traditional clover action is to 
think that the rotations Eq. || are done on the 
lattice JjJ . In this viewpoint the rotations (Eq. ||) 
are similar to the effect of a preconditioner, ex- 
cept that the connection between the twisted and 
standard clover actions is not exact, because the 
transformation in Eq. || also rotates the Wilson 
term in the operator, 

The formalism for the 0{a) improvement of 
the twisted fermion operator has been described 
in§. 

3. INTERPOLATING OPERATORS 

The parity and flavour symmetry violating 
parts of the twisted clover action make the con- 
struction of the interpolating operators nontriv- 
ial. ALPHA g use the twisting rotations in Eq. H 
to construct the interpolating operators. 

Consider a general meson interpolating opera- 
tor 

(5) 

For the Wilson/clover action all three meson op- 
erators (i = 1,2,3) are equivalent. The flavour 
symmetry breaking term in the twisted action 
breaks the equivalence g. 

For i =1,2 the generic meson operators are 

vlTd! = u(cos(a/2) + ?75sin(a/2)) 

r (cos(|) -i7 5 sin(|))d (6) 

For r = 1 (scalar), 75 (pion), 7^ (fei) and 74 7, 
(second rho operator). 

v/Td' = uTd (7) 
For T = 7475 (pion), 74 ("exotic"), 7; (rho) 
u'Td = cos(a)urc? + i sin(a)uj5Td (8) 
The neutral meson operator is proportional to 
TiTu - dTd (9) 



The neutral meson operator has the opposite con- 
nection between the gamma matrix and the mix- 
ing as for the charged meson operators. Discon- 
nected fermion loops are required to compute the 
correlators for the neutral meson operator, be- 
cause the twisted fermion operator does not obey 
Mt= 75 M7 5 [|. 

To compute the correlator for the nucleon with 
an interpolating operator, such as: 

P(x)i = eQbc ( U a (x) T C7 5 d fc (a;)K(x), (10) 

separate inversions are required for the up and 
down twisted clover propagators. These can be 
rotated back to the clover fermion matrix us- 
ing Eq. before the nucleon correlator is con- 
structed. 

4. NUMERICAL RESULTS 

The simulations were performed at = 6.0 
on I6 3 48 lattice in the quenched approximation. 
The numerical value of the clover coefficient de- 
termined by the ALPHA collaboration J^] was 
used. The UKQCD collaboration Q found, at the 
same parameters, evidence for exceptional config- 
urations at Mps/My <~ 0.54, hence we aimed to 
use the twisted clover action to explore the quark 
mass region under that limit. For this exploratory 
study we only used 30 configurations. 

UKQCD found g K cH tical = 0.135252, so we 
initially ran at (k = 0.135, /i = 0.02), (n = 0.135, 
fi = 0.01), and (k = 0.135, fi = 0.005). We were 
aiming for a ~ t|, but this k value was not close 
enough to K cr iti ca i , and the a values for the runs // 
= 0.02, 0.01, and 0.005 were: 0.41tt, 0.33tt, 0.24tt. 
For all these runs we used the standard stabilised 
biconjugate gradient algorithm. To increase a to 
be close to | we ran at (k = 0.13525 , \i — 0.01) 
and [k = 0.13525 , fi = 0.005). For the inver- 
sions with k — 0.13525, the conjugate gradient 
algorithm had to be used. 

With the low statistics we could only obtain 
satisfactory fits to the pion correlator. Factoris- 
ing fits Q were performed to a 2 by 2 matrix, 
with a basis of fuzzed and local correlators ||. 
The pion mass as a function of the renormalised 
quark mass (Mq in Eq. [| with [i q jZp and 
Tdq = ZAfUpcac/Zp, where rripcac is the mass from 
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the PCAC relation ) is plotted in Fig. |[ The pion 
mass from the twisted clover simulation with the 
parameters k = 0.13525,^ = 0.005 is ~ 60 % 
lower than the pion mass obtained from the clover 
action result at the point where exceptional be- 
come a problem (k — 0.13455). 



+ k: Q.'.'.in2ii,/i 0.00b 
« k=0. 13525, ,u=0.01 
O k=0. 13455, ii=0.0 
0.06 he k=0. 135,^=0,005 
□ k = 0. 135, ,u=0.01 

« K=0.135,yU=0.02 



111 of the i^f5/i term on the spec- 
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(12) 



perturbation 
trum is 

d\{p) 
d/i 

where y is the left eigenvector of the clover 
fermion operator Af . Unfortunately the symme- 
tries of the clover operator M do not constrain the 
matrix element of y^jsx. The complex eigenval- 
ues of Af are changed at 0(/i) by the 1/^75 term. 

We are using the ARPACK package to numer- 
ically study the eigenvalues of the twisted clover 
action. 

6. CONCLUSIONS 

In the quenched approximation, the twisted 
clover action has allowed us to simulate at lighter 
pion masses than could be reached by the clover 
action. We are increasing the statistics at the 
lightest quark mass. 

We thank S. Sint and R. Frezzotti for discus- 



Figure 1. The pion mass squared as a function of 
renormalised quark mass. 

5. ANALYSIS OF THE EIGENVALUES 

It is instructive to consider the effect of the 
75 term on the eigenvalue spectrum of the clover 
fermion matrix. Consider the eigenvalue equation 
for a single flavour of the twisted clover action. 



(M + ifij 5 )x(p) = \{/j,)x(h) 



(11) 



where Af is the fermion matrix of the clover ac- 
tion. The relation Aft = 75M75 implies that the 
eigenvalues of Af are either real, or come in com- 
plex conjugate pairs (see Jl(J for a review). 

The real eigenvalues of Af are also eigenvalues 
of 75 , hence the eigenvectors of the clover fermion 
matrix with real eigenvalues are also eigenvectors 
of the twisted matrix with complex eigenvalues. 

The complex eigenvalues of Af are not eigenval- 
ues of 75 , so there is no simple connection between 
the eigenvalues of the twisted clover fermion op- 
erator and those of the clover fermion operator. 
Perturbation theory shows that the first order 
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